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Abstract 

We investigate the power-suppressed corrections to the mean values of various 
quantities that characterise the shapes of final states in deep inelastic lepton 
scattering. Our method is based on an analysis of one-loop Feynman graphs 
containing a massive gluon, which is equivalent to the evaluation of leading 
infrared renormalon contributions. As in e + e~ annihilation, we find that the 
leading corrections are proportional to l/Q. We give quantitative estimates 
based on the hypothesis of a universal low-energy effective coupling. 
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1 Introduction 



The study of final-state properties in deep inelastic lepton scattering (DIS) has received a 
great impetus from the increasing quantity and kinematic range of the HERA data. The 
determination of the strong coupling a s from final-state properties in DIS is an attractive 
possibility because of the relative simplicity of the lepton-hadron interaction, combined 
with the wide range of dynamical scales available in a single experiment at a single beam 
energy. It is expected that the scale for a s will be set primarily by the lepton-hadron 
momentum transfer-squared Q 2 , which can range from zero to 10 5 GeV 2 at HERA. Thus 
there should be a wide region in which the value and running of a s (Q 2 ) can be observed 
with good precision. 

One possible method for a s determination is the measurement of jet fractions Jl|, defined 
according to one of the several available infrared-safe jet algorithms By definition, 

jet rates defined by an infrared-safe algorithm can be computed in perturbation theory, 
and next-to- leading-order calculations are now available ||||. The a s values obtained by 
comparing jet rates with HERA data are consistent with those found in other processes, 
and in particular they show the expected decrease with increasing Q 2 . 

In the present paper we consider a different set of DIS final-state observables which can 
be used to determine a s . These are the various event shape variables which can be defined in 
analogy with those used in the study of e + e~ annihilation final states. In e + e~ physics, event 
shapes have been found to be a useful tool for testing QCD and measuring a s . They can 
be defined so as to be sensitive to different aspects of QCD dynamics (e.g. the longitudinal 
or transverse development of jets) are subject to different non-perturbative 'hadronization' 
corrections. Thus a s determinations from a variety of event shapes complement those from 
jet rates and give an indication of the systematic uncertainties due to non-perturbative 
effects. The same considerations make it important to calculate and measure event shapes 
in DIS. 

Another reason to measure event shapes is that there are new theoretical ideas about 
non-perturbative corrections to them []6| [10|1 , which provide constraints on a s determina- 
tions from event shapes and are interesting to test in their own right. By looking at the 
behaviour of the QCD perturbation series in high orders, one can identify unsummable, 
factorially divergent sets of contributions (infrared renormalons JTTJ] ) which indicate that 



non-perturbative power-suppressed corrections must be included. The Q 2 -dependence of 
the leading correction to a given quantity can be inferred, and by making further uni- 
versality assumptions one may also estimate its magnitude. Tests of these ideas provide 
information on the transition from the perturbative to the non-perturbative regime in QCD. 
In particular, one can investigate the possibility that an approximately universal low-energy 
effective coupling may be a useful phenomenological concept [|T|,|5|,[r2] . 

Such an approach has been applied with some success to e + e~ event shapes [f7i[8|Jl3i 
and fragmentation functions WM, and to DIS structure functions p|JT3|. In the present 



paper we extend it to event shapes in DIS JT6|. We find that, as in e + e annihilation 



the leading power corrections to these quantities are typically proportional to 1/Q. The 
hypothesis that they are related to a universal low-energy effective coupling implies that 
their magnitudes are given by a single non-perturbative parameter. We give quantitative 
estimates based on the value of this parameter derived from e + e~ data. 
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In the following Section we explain how the DIS event shape variables that we com- 
pute are defined. In Sect. 3 we give the leading-order perturbative predictions for these 
quantities. To determine a s , one needs the predictions in next-to-leading order, which are 
not yet available. However, the leading-order calculation provides a useful guide to the 
relative importance of the power-suppressed corrections, which we estimate in Sect. 4 us- 
ing the method of Ref. f|. We explain how these estimates can be refined and combined 
with the next-to-leading predictions when they become available. Finally, our results are 
summarized briefly in Sect. 5. 



2 Event shape variables in DIS 

A complication in DIS, absent from e + e~ annihilation, is the presence in the final state 
of the remnant of the initial-state hadron, i.e. the constituents that did not participate in 
the hard scattering of the lepton. It is expected that the fragmentation of the remnant 
will be dominated by soft, non-perturbative physics. While of interest for studying the 
hadronization process, the remnant fragmentation is not so useful for a s determinations, 
and therefore we concentrate here on aspects of event shapes that are not sensitive to 
it. This is conveniently done by looking at the final state in the Breit frame of reference 

OSS- 



We consider the deep inelastic scattering of a lepton of momentum / from a nucleon of 
momentum P, with momentum transfer q. The main kinematic variables are Q 2 = —q 2 , 
the Bjorken variable x = Q 2 /2P ■ q and y = P ■ q/ P ■ I ~ Q 2 /xs, s being the total cm. 
energy squared. Then the Breit frame is the rest-frame of 2xP + q. In this frame the 
momentum transfer q is purely spacelike, and we choose to align it along the +z axis: 

P = lQ(l/x, 0, 0, -l/x) , q = ±Q(0, 0, 0, 2) . (2.1) 

To a good approximation, the fragmentation products of the remnant will be moving 
in directions close to that of the incoming nucleon, i.e. they will remain in the 'remnant 
hemisphere' H r (p z < 0). On the other hand the products of the hard lepton scattering will 
tend to be found in the 'current hemisphere' H c (p z > 0). In fact in the parton model the 
scattered parton moves along the current (+z) axis with momentum xP + q = |Q(1, 0, 0, 1). 
Thus in the parton model the current hemisphere looks like one hemisphere of the final 
state in e + e~ annihilation at centre-of-mass energy Q. Fragmentation studies have shown 
that this similarity is indeed manifest in hadron spectra and multiplicities ||20|| . This makes 



it natural to define event shape variables in close analogy to those for e + e annihilation, 
but limited to particles a appearing in the current hemisphere, a G H c . 

We can now construct infrared-safe quantities that characterize the shape of the event 
defined in this way. Perhaps the simplest is the current jet thrust [[HJ 

T Q = 2 £ p a ■ n/Q (2.2) 

aeH c 

where n represents the unit 3-vector along the current direction (the +z axis, in our 
convention). The subscript Q indicates that T is normalized to \Q. Alternatively we may 
normalize to the total energy in the current hemisphere, 

Te = E Pa ■ n/ £ E a . (2.3) 

a£H c a£H c 
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Both of these quantities are equal to unity in the Born approximation, and their deviation 
from this value measures the longitudinal development of the current jet. It will therefore 
be convenient to study instead the quantities tq = 1 — Tq and t e = 1 — T E , which vanish 
in the Born approximation. 

It is kinematically possible for the Breit frame current hemisphere to be empty. In that 
case, taken literally, Eq. (|2.2|) implies that Tq = 0, hence tq — 1, while Eq. ( |2.3| ) leaves r E 
undefined. For consistency with the other event shapes defined below, we instead define 
t q — t e = when the current hemisphere is empty. 

Similarly we can define the current jet broadening [[H 



Bq=J2 \p a xn\/Q, (2.4) 



aeH c 

or 



B E = \Y.\Pa* n\/ £ E a , (2.5) 

aeH c aeH c 

which emphasizes the transverse development of the jet. 

Both the thrust and the broadening are defined here with respect to the current direction 
n. Two quantities which measure the jet development independent of direction (apart from 
the restriction to particles in the current hemisphere) are the scaled current jet mass 



(2.6) 




and the C -parameter | 22fl 



C Q = 3(AiA 2 + A 2 A 3 + A 3 Ai) (2.7) 
where Ai^.3 are the eigenvalues of the linearized momentum tensor 



0« = 2 £ (pX/|pj)/Q. (2.8) 



aeH c 

Again, we may alternatively define quantities p E and C E , in which Q is replaced by twice 
the total energy in the current hemisphere, so that 

Pe/pq = C E /C Q = Q 2 / Uj2E a } . (2.9) 

\ aeH c 



3 Leading-order perturbation theory 



At first order in a s , up to two final-state partons can be emitted in the hard lepton-parton 
subprocess, as illustrated in Fig. 0. The momentum of the struck parton is p = xP/^ 
(x < £ < 1) and we define z = P ■ r/P ■ q (0 < z < 1). 

The differential cross section is 



d 3 a 



2vra 2 



dxdQ 2 dz Q 4 



{ [l + (1 - yf] F T (x, z) + 2(1 - y)F L (x, z)} 



(3.1) 
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Figure 1: Jet production in deep inelastic scattering. 



The generalized transverse and longitudinal structure functions Ft(x,z) = 2F\(x,z) and 
F L (x, z) = F 2 (x, z)/x — 2F 1 (x, z) are of the form (for z < 1) 



Fi(x, z) = ^ ^ [CpC^Ce, z)q(x/0 + TfC iig (^, z)g(x/0] 



where 



q(x) = $> 2 [Qj( x ) + Qj( 

3=1 



X 



3=1 



for / active quark flavours, Cp = 4/3, Tr = 1/2 and JT8 

£ 2 + z 2 



^'^-(i-0(i-,) 

C x , ? (£,z) = 4£z 

Pr^«) = [e+(i-o 2 \ 

C L j{Z,z) = 8Z(l-Z). 



+ 2£z + 2 



z 2 + (l- z) 

2(1 - *) 



In the Breit frame P and g are given by Eq. (|2.1| ) and we can write 

P = ig(i/e,o, o,-i/o 

r = |Q(^o, 0, z 3 ) 



where 



fc = \Q{zq, -zx,0,z 3 ) 
z = 2z-l + (l-z)/£ 

z 3 = 1 - (1 - ^)/e 

2 = 1 - 2^ + z/f 

^3 = 1" */£ 



z ± = 2^(1 (!-£)/£■ 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 
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Figure 2: Phase space region for jet production in deep inelastic scattering. 



We can distinguish four subregions of phase space, as illustrated in Fig. ^: 

A: both produced parton momenta k, r in the current hemisphere (z^,z 3 > 0); 

B: only parton momentum r in the current hemisphere (2:3 > 0, Z3 < 0); 

C: only parton momentum k in the current hemisphere (z 3 < 0, Z3 > 0); 

D: no produced parton momenta in the current hemisphere (2:3, z 3 < 0). 

In leading order the event shape variables defined in Sect. are given in these regions 
by Table 1. By definition they are all zero in region D. By construction, they all vanish in 
the soft and/or collinear limits £, z — > 1. Note that p and C also vanish throughout regions 
B and C. 

At any particular values of x and Q 2 , the mean value of a shape variable S is now given 
in leading order by 



where 



and the denominator is the differential cross section evaluated in Born approximation, 

d 2 a 2ira 2 



dxdQ 2 " Q 4 



l + (l-yy q(x). (3.9) 



We discuss the numerical values of the leading-order predictions together with the power 
corrections in the following Section. 
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Table 1: Event shape variables S(£, z) in leading order. 
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4 Power corrections 



Our estimate of the leading power corrections to the perturbative results given above is 
based on the approach of Ref . || . Non-perturbative effects at long distances are assumed 
to give rise to a modification Sa e s(fi 2 ) in the QCD effective coupling at low values of the 
scale /i 2 . The effect on some observable F is then given by a characteristic function T{x, e), 
as follows: 

da 2 

J^5a e{i (fi 2 )F(x,e = v 2 /Q 2 ) (4.1) 



SF(x,Q 2 



/r 



where 



F{x,e) 



(4.2) 



The characteristic function is obtained by computing the relevant one-loop graphs with a 
non-zero gluon mass /i = Q-Je |||23]. 



Arbitrary finite modifications of the effective coupling at low scales would generally 
introduce power corrections of the form 1 / /x 2p into the ultraviolet behaviour of the running 
coupling a s itself. Such a modification would destroy the basis of the operator product 
expansion fl24 |. One must therefore require that at least the first few integer moments of 
the coupling modification should vanish: 

dfj, 2 



n c 



{j/f Sa cS (fi 2 



0; p=l, 



,Pn 



(4.3) 



The upper bound p max could be set by instanton-anti-instanton contributions (p max ~ 9). 
The constraint ( |4.3j ) means that only those terms in the small-e behaviour of the character- 
istic function that are non- analytic at e = will lead to power-behaved non-perturbative 
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contributions. These are just the the terms that give rise to infrared renormalons in per- 
turbation theory |23| . 



For gluon mass-squared /i 2 = eQ 2 the quark coefficient functions in Eq. ( |3.4|) become 
r ,t , _ (1 ~ z)(l - Q + 2fr(l - zf - ft 2jz{\ - e) (1 - z)(l - Q - £e 

4^(1 - z) 2 
^'^(l-.-ft) 2 ' 

(4.4) 

Processes involving an incoming gluon are not expected to give terms that are non-analytic 
at e = 0, and therefore we do not consider them as a source of power corrections. The 
kinematic variables that give the momenta according to Eq. Q3.5|) are now 



Zq = 2z-1 + (1- z)/£ - e 
z 3 = l-(l-z)/£ + e 

z = 1 - 2z + z/i + e (4.5) 
z 3 = 1 - zji- e 



z x = 2 s Jz{l-z){l-i)/i-ez. 

Thus the phase space region is now < z < 1 — e£/(l — £), as illustrated in Fig. |3]. The 
regions A,. . . D defined above in terms of the signs of Z3 and z% are as indicated. 

The corresponding characteristic function for the mean value of some event shape vari- 
able S is given by Eq. (|3.7|) with F^ S \x) replaced by (Cf/2 / k)J^ S ^(x, e) where 

rt S \x,e)= f 1 ^ f 1 dzS^,z,e)C iiq (i,z,e)e(l-z-^ + ^z-eOq(x/0- (4.6) 

Note that for brevity we have extracted the overall factor of Cp/2ir. The expressions for 
the shape variables S(£, z, e) are as given in Table 1, but with the kinematic variables now 
given by Eqs. (|4.5|) instead of Eqs. 
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Finding the leading non- analytic term in the behaviour of the integral ( f4.6|) as e — > 0, 
differentiating as instructed in Eq. ( }4.2| ), and inserting the result in Eq. ( f4.f| ), we obtain 
the corresponding predicted power correction. Generally speaking, the small-e behaviour 
of the characteristic function is dominated by the region around the boundary point P in 
Fig. ^], and therefore the leading power correction is independent of whether we normalize 
the shape variable to Q/2 or to the energy in the current hemisphere. 

Following Ref. ||, we may express the magnitudes of power corrections in terms of the 
moment integrals 

A 2P = ^ r^-fj^Sa^), (4.7) 

27T JO fa 

which vanish for integer p, and their p-derivatives 



=^r% ^ ^ <w/i 2 ) , (4.* 

2,71 Jo fa 



which are in general non- vanishing for any p. The leading corrections to the event shapes 
we are considering correspond to p — ~, and therefore, on the assumption that 5ot e R is 
universal, they can all be expressed in terms of the two non-perturbative parameters A\ 
and A[. Studies of event shapes in e + e~ annihilation suggest that A\ ~ 0.25 GeV ||, with 
A[ as yet undetermined. 

As an alternative representation of the magnitudes of power corrections, we may adopt 
the approach of Ref. J!]] and express them directly in terms of moments of a e g over the 
infrared region. We substitute for 5a e g in Eq. (|4.7|) 

5a cS (/2 2 ) ~ a cfT (/i 2 ) - ttfV) , (4.9) 

where af T represents the expression for a s corresponding to the part already included in 
the perturbative prediction. As discussed in Ref. f7j, if the perturbative calculation is 
carried out to second order in the MS renormalization scheme, with renormalization scale 
fa\, then we have 

afV) = « s (/4) + ibM^i/v 2 ) + k] « s 2 (/4) (4.10) 

where (Ca = 3) 

llCi-2/ (67-^)0,-10/ 

127T 367T ^ J 

The constant k comes from a change of scheme from MS to the more physical scheme []25| in 
which a e ff is defined. Then above some infrared matching scale fa we assume that a e fj(/i 2 ) 
and af T (/x 2 ) approximately coincide, so that 



M - TT I 1 («cfr(/i 2 ) - a s (/i 2 ) - [Mn(/4/V) + k] a 2 s {fa 

2tt Jo fa v 



F ^ 1 d ^j-V> (a eS (v 2 ) -a s (nl) - \bMnllij 2 ) + k]a 2 (^ 2 ' 

C 

— fa (a (fa) - ot 8 (fj.l) - [Mn(/z 2 //i 2 ) + k + 2b] a 2 s (fa\] 

7T x 



(4.12) 



where 



1 f« „ 

a (fa) = — / oi eS (n ) d\x . (4.13) 
fa Jo 

Thus in this notation the value of Ai determines the average value of the effective coupling 
below the matching scale fa. The dependence of ao on fa is partially compensated by the 
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/^-dependence of the other terms on the right-hand side of Eq. ( 4.12| ). The dependence 
on the renormalization scale should help to compensate the scale dependence of the 
perturbative part. Notice that if we take /i| oc Q 2 then Ai has a logarithmic dependence 
on Q 2 . In general we do expect 'power' corrections to have additional logarithmic Q 2 - 
dependence (anomalous dimensions), but this cannot yet be calculated reliably for event 
shapes. 

In Ref. it was found that the formula (|4.12|) with /4 = Q 2 , /i : = 2 GeV and a (2 GeV) 
= 0.52 gave good agreement with e + e~ event shape data. Similar results were obtained in 



Ref. [13 



4.1 Current jet thrust 

In the case of the shape variables tq or tb, the behaviour of the expression (|4.6|) for the 
transverse contribution Tt\ 

e — > is found to be of the form 

^\x,e) ~ f£\x,0) -%yf~eq{x) , (4.14) 

while the longitudinal part J 7 ^ is less singular at e = 0. Thus from Eqs. ( |3.7| ), (|Q| ) and 
( |4.7|) we obtain the leading non-perturbative contribution 

*<t>~4^. (4.15) 

The behaviour (|4.14| ) at small e follows from the fact that the derivative is dominated 
by the phase space boundary z = l — e£/(l — £): 

4 T) (a;, e)~ £ f ^ /' dz r(£)C^(£, z, 0)6(1 - z - £ + £z - e£) q(x/0 . (4.16) 

Jx JO 

From Table 1 and Eq. (|4.5|), on this boundary we have 



r(e,z,0) = (l-O/^ fore>^p 
= (!"*)/£ for^^P 



(4.17) 



where = 1/(1 + y/e). Thus 



. (T) rtpdt i+e (x\ r^d^i+e fx\ r 

in agreement with Eq. Q4.14Q . 

Numerical predictions for the mean value of the current jet thrust in ep scattering at 
y/s = 296 GeV are shown in Fig. [| as a function of Q for various values of x. The MRS 
A' parton distributions [26] were used, with the corresponding value =231 MeV in 
the two-loop expression for a s {Q 2 ). The leading-order perturbative predictions given by 
Eq. ( |3.7| ) are shown by the dashed curves. For the power correction coefficient A% we 
used Eq. Q4.12j ) with /i| = Q 2 , Hi = 2 GeV and a (2 GeV) = 0.52, as in the fits to e+e" 
data, but we omitted the term of order a 2 because we are combining with only a first-order 
perturbative calculation in this paper. When the higher-order prediction becomes available, 
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Figure 4: Predictions for the mean value of the current jet thrust in deep inelastic scattering. 
Left- and right-hand plots are for the two definitions (2.2) and (2.3), respectively. Dashed: 
leading-order perturbation theory. Solid: leading order plus leading power correction. In 
each case the four curves (top to bottom) are for x = 0.003, 0.01, 0.03, 0.10. 



the 0{a 2 g) term in Eq. (|4.12| ) should be included when estimating the power correction to 
it. 

The resulting overall predictions are shown by the solid curves. We see that the esti- 
mated power correction is substantial, 20-30% at Q = 50 GeV and dominating below 15 
GeV. There is significant x-dependence in the perturbative prediction, and large- X data 
(x > 0.01) are required to cover the region where the power correction is under control. 
With sufficient data in the range Q — 15 — 50 GeV, however, it should be possible to 
perform a two-parameter fit to determine a s and cxq from the average current jet thrust. 



4.2 Current jet broadening 

For the jet broadening Bq or Be we find a slightly different behaviour at small e, namely 

^ B) (x,e) ~ j4 B) (x,0) + 8v^(lne + c)g(x) , (4.19) 

where c is a constant (probably x-independent) which we cannot determine reliably. This 
implies a non-perturbative correction of the form 

<5(B)~4^(lnQ 2 -c-2)-4^-. (4.20) 

Since we do not know the value of c, we may as well absorb all the non-logarithmic terms 
into an unknown scale, Qq\ 

5(B) = 8^\n(Q/Q ). (4.21) 
10 
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Figure 5: Predictions for the mean value of the current jet broadening in deep inelastic 
scattering. Left- and right-hand plots are for the two definitions (2.4) and (2.5), respectively. 
Curves as in Fig. 3. 



We find that, in contrast to the situation for the thrust, the result ( |4.19| ) is only obtained 
when one includes the gluon mass explicitly in the definition of the jet broadening. This is 
because, unlike the thrust, the broadening vanishes on the phase-space boundary z± = 0. If 
we neglect the gluon mass in the definition, this is not the case and an expression analogous 
to Eq. (|4.16|) is obtained, which gives 



We then find 



j^\e) ~ -3^(lne + c) , (4.23) 

corresponding to a coefficient of 6 instead of 8 in Eq. ( [4.211) . Thus the correct form is 
obtained, but the full mass-dependence must be retained to compute the coefficient. 

The fact that the magnitude of the leading power correction to the jet broadening is 
sensitive to the gluon mass- dependence in its definition suggests to us that the prediction 
for this shape variable less reliable than that for the thrust. As pointed out in Ref. |J, 
shape variables are not fully inclusive with respect to the fragmentation of the gluon: their 
values for the 'decay products' of a timelike virtual gluon are not necessarily equal to those 
for a 'real' gluon of equivalent mass. In the case of the thrust, model studies suggest that 
the numerical effect of this on the leading power correction is small, but we expect it to be 
larger for variables that depend explicitly on the gluon mass. 

Numerical predictions for the current jet broadening at HERA are shown in Fig. |j| using 
the same parameter values as before to compute A\ and, for definiteness, Qq = /x r = 2 GeV 
in Eq. ( 4.21Q . We see that the resulting power corrections are large. As we have stressed 
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Figure 6: Predictions for the mean value of the current jet mass in deep inelastic scattering. 
Left- and right-hand plots are for the two definitions (2.6) and (2.9), respectively. Curves 
as in Fig. 3. 



above, they are also more uncertain in this case, suggesting that jet broadening in DIS is 
not a good shape variable for a s determinations. 



4.3 Current jet mass 

Next we consider the power corrections to the jet mass pq or p#. We notice from Table 
1 that there is no explicit gluon mass dependence in the definition of these variables, and 
there is no contribution outside the phase space region A.Q Inside this region we have in 
fact pq = tq. Thus the jet mass receives a contribution from the second integral only in 



Eq. ( 4. IS ). This gives exactly one half of the correction to the thrust and so one finds that 

J#\x, e) ~ J$\x, 0) - A^q(x) (4.24) 
at small e, which implies a non-perturbative correction 

5(p)~2^. (4.25) 



The numerical predictions for the current jet mass, shown in Fig. |6], suggest that this 
is a good variable for a s determinations. The power correction is somewhat larger than 
that for the thrust, relative to the perturbative prediction (cf. Fig. but there is less x 
dependence in the latter. 

' There is a contribution e in region C, but since this is analytic it does not contribute to the power 
correction. 
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Figure 7: Predictions for the mean value of the C-parameter in deep inelastic scattering. 
Left- and right-hand plots are for the two definitions (2.7) and (2.9), respectively. Curves 
as in Fig. 3. 



4.4 C-parameter 

Finally we compute the power correction to Cq or Cg. Here again we see from Table 1 
that that there is no contribution outside the phase space region A. However in this case 
the variable, unlike the jet mass, does depend explicitly on the gluon mass. From a full 
evaluation retaining this mass dependence we find the small-e behaviour 

T ( t\x, e) ~ FP(x, 0) - 2An^q(x) , (4.26) 

corresponding to a leading non-perturbative contribution 

5 (C) ~ 12tt^- . (4.27) 

If one uses the massless gluon expression for C, one obtains 

T>n ' i, e) ~ 12e / — ; — ; ~ 6nWe , 4.28 

TK } A/(i+Vi) e l-f + e(l-20 (l-0 2 + ee(2^-l) V ' 1 ' 

which corresponds to one-half of the full result. We would argue again that this sensitivity 
to the gluon mass- dependence of the definition suggests that the prediction for the C- 
parameter is less reliable than that for the thrust and jet mass. 

The numerical results, Fig. [7], show that the estimated power correction is very large in 
this case, dominating over the perturbative prediction even at Q = 50 GeV. The size and 
uncertainty in the correction suggest that, like the jet broadening, the C-parameter is not 
a good variable for determining a s . 
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5 Summary 



In this paper we have investigated several infrared- safe variables which characterize the 
shapes of DIS final states in the current hemisphere of the Breit frame, where one avoids 
as far as possible complications associated with the target remnant. We have presented 
numerical predictions to leading order in perturbation theory, together with estimates of 
leading non-perturbative power corrections, which are predicted to be proportional to 1/Q, 
modulo logarithmic Q-dependence. The assumption of an approximately universal low- 
energy effective coupling allowed us to relate the magnitudes of the corrections to those 
in e + e~ annihilation. We found that they are expected to be largest, and most uncertain, 
for the current jet broadening and C-parameter, and so these observables are probably 
not suitable for determination of the perturbative strong coupling a s . The current jet 
thrust and mass should have power corrections that are smaller and under better control. 
When higher-order predictions for these quantities are available, our predictions of the 
power corrections can also be refined, and it should be possible to measure both ot s and 
the relevant non-perturbative parameter, «o(/ii) in Eq. (|4.12| ), from these quantities. 
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